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Overview of Genomic Insights into Chicken Growth Traits Based on Ge- 
nome-Wide Association Study and microRNA Regulation 
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Abstract: Over the two past decades, a significant number of studies have observed animal growth traits to examine ani- 
mal genetic mechanisms due to their ease of measurement and high heritability. Chicken which has a significant impact 
on fundamental biology is a major source of protein worldwide, making it an ideal model for examining animal growth 
trait development. The genetic mechanisms of chicken growth traits have been studied using quantitative trait loci map- 
ping through genome-scan and candidate gene approaches, genome-wide association studies (GWAS), comparative ge- 
nomic strategies, microRNA (miRNA) regulation of growth development analysis, and epigenomic analysis. This review 
focuses on chicken GWAS and miRNA regulation of growth traits. Several recently published GWAS reports showed 
that most genome-wide significant single nucleotide polymorphisms are located on chromosomes 1 and 4 in chickens. 
Chicken growth, particularly skeletal muscle growth and development, is greatly regulated by miRNA. Using dwarf and 
normal chickens, let-7b was found to be involved in determining chicken dwarf phenotypes by regulating growth hor- 
mone receptor gene expression. 
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1. INTRODUCTION 

The growth and development of livestock is a compli- 
cated life process and is subject to a significant amount of 
genetic control. Several studies have examined "hotspot" 
areas, including candidate gene mutation, effects of quantita- 
tive trait loci (QTL), genome-wide association studies 
(GWAS), microRNA (miRNA) regulation, and epigenetic 
inheritance. Gene expression patterns can be used to explain 
phenotype variations and to illustrate epigenetic mechanisms 
such as consensus methylation-rich promoters that corre- 
spond to down-regulation of gene expression and transgen- 
erational stability of epigenetic variants that simultaneously 
appeared in mammals, birds, and plants [1-3]. Focusing on 
growth traits and miRNA regulation, this review discusses 
three main areas related to chicken growth. 

First, growth trait-correlated mutations in the genes of the 
somatotropic axis (Table 1) will be discussed before focus- 
ing on GWAS and miRNA regulation of growth traits. A 
number of single nucleotide polymorphisms (SNPs) located 
within genes encoding for growth hormone (GH) and its 
receptor (GHR), ghrelin (GHRL) and its receptor (GHSR), 
insulin-like growth factor 1 (IGF-1) and its receptor (IGF- 
1R), insulin (INS), insulin-like growth factor binding protein 
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2 (IGFBP2), adipose triglyceride lipase (ATGL), and pitui- 
tary-specific transcription factor 1 (PIT1) have been associ- 
ated with chicken growth traits [4-17]. Moreover, SNPs of 
the high mobility group AT-hook 2 gene (HMGA2), which is 
located on chicken (Gallus gallus) chromosome (GGA) 1, 
should be noted for their significant associations with body 
weight. Common insertion-deletion (indel) variations in 
growth-related genes have also been shown to act on growth 
diversity. The 1,773 bp deletion located at exon 10 and the 
3' untranslated region (UTR) of the GHR gene resulted in a 
dysfunctional precursor, leading to sex-linked dwarf chick- 
ens (SLD). One 57-bp indel in intron 2 of the PIT1 gene was 
reported to be significantly associated with hatch weight and 
shank length at 84 d of age [8]. Another 6-bp indel of the 
GHSR gene was found to be related to crude fatty content of 
the leg muscle, although no evidence for its association with 
growth traits was identified [7]. An 8-bp indel in exonl of 
the GHRL gene may have a negative effect on chicken 
growth [13]. Some SNPs are likely fixed and must transmit 
in a high frequency format through haplotypes such as strong 
linkage disequilibrium. Using haplotypes and linkage dise- 
quilibrium, important SNPs and QTL were identified. For 
instance, no indication of an association between either SNP 
A17299834G or SNP C17293932T and the IGF-IR gene 
was observed for growth traits; however, haplotypes could 
be constructed from these two SNPs, which were shown to 
be significantly associated with body weights, daily weight 
gains, and leg length [9]. Additionally, SNP C6540334T and 
SNP C654201 IT of the GHR [10] gene and SNP A+428G of 
the INS gene [14] exhibit a similar pattern. 
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Table 1. Significant Associations of the SNPs Within Candidate Genes with their Growth Traits (P<0.01) 
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Gene 


Chr 


SNP" 


Growth Traits b 


Ref. 


GHRL 


7 


C71T 

G1215A 
8bp indel (exon) 


BW: 7 wk, 16 wk; 
SG: 7 wk, 16 wk; SL: 16 wk 
BW: 16 wk 

BW: 2 wk, 3 wk, 5 wk, 6 wk, 7 wk, 13 wk; BL 


11 

13 


INS 


5 


T+3737C (intron) 
A+3971G(intron) 
C+1549T (intron) 


SIL 

BW: 4 wk; HW; BA; SIL 
BW: 4 wk, 12 wk; BD 


14 


GH 


1 


G+119A (intron) 
G+1705A (intron) 
G+3037T (intron) 


SIL 

BW: 3 wk, 4 wk, 10 wk, 12 wk; ADG: 0-4 wk; 
SL: 10 wk, 12 wk 


15 


IGF1R 


10 


A17299834G (5'UTR) 


LL: 8 wk 


9 


GHR 


z 


G6631778A (3'UTR) 


HW; BW: 5 wk, 6 wk, 7 wk 


10 


PIT1 


1 


57 bp indel (intron) 
rsl3687127 (intron) 
rsl3687128 (exon) 


HW; SL: 12 wk 
SD: 1 1 wk 
BW: 4 wk; ADG: 0-4 wk 


8 


ATGL 


5 


c.782G>A (exon) 


BW: 6 wk, 7 wk, 8 wk 


6 


IGFBP2 


7 


G729T (intron) 
C1032T (intron) 
A663T (exon) 
G738A (exon) 


BW: 7 wk 
HW; BW: 2 wk 
BW: 1 wk, 2 wk, 3 wk, 5 wk, 6 wk, 7 wk; HW 
BW:1 wk, 2 wk, 3 wk, 4 wk, 5 wk, 6 wk, 7 wk, 8 wk, 13 wk; HW 


16 


RBI 


1 


g.39692 G>A 
g.77260 A>G 


BW: 6 wk, 1 1 wk 
BW: 5 wk, 6 wk, 7 wk, 8 wk, 9 wk, 10 wk, 1 1 wk, 12 wk 


24 


HMGA2 


1 


rs 1523 1472 (intron) 
rs 13 8493 81 (intron) 


BW:6 wk 

BW: 1 wk, 2 wk, 3 wk, 5 wk, 6 wk; HW 


5 



8 Names of SNPs are cited from original papers. 

h BW: body weight; ADG: average daily gain; SL: shank length; BL: body length; LL: leg length; SD: shank diameter; HW: hatch weight; SIL: length of small intestine; BA: breast 
angle; BD: breast depth; SG: shank girth. 



In addition to candidate gene approaches, identifying QTL 
underlying traits is a useful strategy. Continuous improvement 
in molecular genetics technology and DNA marker availabil- 
ity, identifying QTL controlling growth traits for marker- 
assisted selection applications has rapidly increased [18, 19]. 
More QTL have been identified on GGA1 and GGA4 than on 
other chromosomes. Among growth traits correlated with 
known QTL, body weight accounts for 30% (207), growth rate 
for 6% (40), and abdominal fat weight for 5% (35), while 
other traits constitute the remaining traits (http://www. animal 
genome.org/cgi-bin/QTL db/GG/index). For GGA4, most 
reported QTL are growth-related. Although several QTL have 
been successfully identified, most are orientated in dispersive 
regions of distinct chromosomes. They typically span a large 
range. For example, QTL affects abdominal fat weight on 
GGA1, which ranges from 73.6 to 215.4 cM [20] and 454 to 
276 cM [21]. This makes applying these results to various 
commercial lines infeasible because imprecise mapping. In- 
creasing marker density or population size under tightly re- 
stricted conditions is used for trouble-shooting. In a study by 
Liu et ah, 10 QTL were identified at the 1% chromosome- 
wide level for body weight from 4 to 12 weeks [22]. Shortly 
after, the authors refined the previously reported QTL interval 
by increasing nine novel microsatellites between LEI0079 and 



ROS0025 and enlarging family sizes from 4 to 12 and F2 in- 
dividual numbers from 369 to 1011, leading to more precise 
QTL identification for these two traits across this region [23]. 
For fine mapping of this QTL, 14 additional SNPs close to 
ADL328 were examined and the most likely QTL for body 
weight was orientated between SNP G173196313A and SNP 
T173539900C over a 400-kb interval [24]. Finally, two SNPs 
(g.39692 G>A and g.77260 A>G) of the retinoblastoma 1 
(RBI) gene were confirmed to be related to growth. Using 
high-resolution QTL is effective for identifying quantitative 
traits near genes. 

2. GWAS ON CHICKEN GROWTH TRAITS 

The genome-wide association study, which shows a re- 
markable success rate, can be used to detect sequence varia- 
tions influencing growth traits. GWAS on chicken growth 
traits showed that most significant SNPs were centered near 
GGA1 and GGA4 [25]. Growth traits including aggregate 
body weight at 0-90 d of age measured weekly, biweekly 
average daily gains derived from weekly body weight, and 
breast muscle weight, leg muscle weight, and wing weight at 
90 d of age were identified in a 1.5-Mb region (173.5-175 
Mb) of chicken GGA1, which was a narrow range compared 
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to other SNPs evaluated. Of effective SNPs, five SNPs in the 
1.5-Mb karyopherin alpha 3 (KPNA3)-forkhead box Ol 
(FOXOIA) regions showed the highest significant effects for 
all growth traits, while two miRNA genes putatively target- 
ing mRNA of IGF-1, FOXOIA, and KPN A3 genes were con- 
tained. Another GWAS study showed that most significant 
SNPs associated with late growth during weeks 7-12 were 
present in an 8.6-Mb region (71.6-80.2 Mb) of chicken 
GGA4, whereas only one SNP on GGA18 was significantly 
associated with early growth (body weight at 2 w of age). 
This may be due to that more genetic variance during early 
growth is caused by epistatic interactions rather than by sin- 
gle point effect [26, 27] . According to the results of artificial 
selection revealed in another study, a region (60-80 Mb on 
GGA4) may have been under intense recent selection in di- 
vergent chicken lines, where 50 generations of selection 
showed 9-fold differences in body weights. Similarly on 
GGA1, several regions respond highly to selection [28]. This 
implies that these two confined regions are important for 
understanding the molecular basis of growth traits. Although 
a lack of conformity between the two GWAS studies was 
caused by breed diversity (different materials), previously 
reported QTL offer powerful evidence supporting that these 
two regions and chicken growth are related (http://www. 
animalgenome.org/cgi-bin/QTLdb/GG/index). Because the 
described DNA variants identified through GWAS may in- 
fluence growth traits, subsequent experiments should be per- 
formed [29, 30]. Of the 15 validated SNPs in the region from 
173.5 to 175Mb, 12 SNPs with P<0.0033 were significant 
between fast-growing and slow-growing groups, indicating 
that these SNPs play a role in a causal mechanism or in cou- 
pling linkage phase with causal mutations [25]. 

The most promising genes harboring growth-related 
SNPs, including the FOXOIA gene, in which two highly 
significant SNPs at 8.9 Kb upstream and 1.9 Kb downstream 
were identified, may contribute to myogenic growth and 
differentiation. FOXOIA expression is notably up-regulated 
in energy-deprived skeletal muscle. Transgenic mice overex- 
pressing Foxola in skeletal muscle weighed less than wild- 
type mice and showed reduced skeletal muscle mass [31]. 
The integrator complex subunit 6 (INTS6) gene, containing 
the most significant SNP effecting aggregate body weight at 
90 d of age, may act as a tumor suppressor gene by inhibit- 
ing growth of prostate cancer by altering the cell cycle pro- 
file and Wnt signaling, resulting in down-regulated expres- 
sion in human multiple prostate cells [32]. The LIM domain- 
binding factor 2 (LDB2)gene, in which SNPs important for 
late growth are located, can bind to several transcription fac- 
tors and drive vascular maturation [33]. 

SNPs in GGA1 and GGA4, which are reliably predicted 
using powerful statistical methods, can be used to understand 
chicken growth. Phenotypic differences should not only be 
ascribed to genetic differences or epigenetic changes as a 
result of sequence mutations, but also to miRNAs which 
affect this intricate phenomenon. 

3. MICRORNA REGULATION OF CHICKEN 
GROWTH 

3.1. MicroRNA Regulation of Biological Processes 

MicroRNAs are key elements in several cell processes, 
including proliferation, differentiation, and apoptosis. They 



are a class of small non-coding RNAs approximately 22 nu- 
cleotides in length and are generally expressed in different 
tissues, regulating expression of specific genes at the post- 
transcriptional level by targeting mRNAs for degradation or 
inhibiting translation [34]. During miRNA formation, 
miRNA genes are initially transcribed as long primary tran- 
scripts (pre-miRNAs) by RNasell or III. Mature miRNAs are 
generated after processing by Dicer, an RNase III endonu- 
clease, and subsequently incorporate into the RNA-induced 
silencing complex (RISC) [35]. Mature miRNA always 
binds to the 3' UTR of target mRNAs, which mainly induces 
translational inhibition and exonucleo lytic mRNA decay; 
targets can also be cleaved endonucleolytically as a result of 
highly complementary base-pairing [36]. Each miRNA has a 
specific structure and formation mechanism. For example, 
subtle sequence variations exist within the gga-miR-1 and 
gga-miR-133 genes, but none exist for gga-miR-206 
(http://www.mirbase.org/). Gga-miR-la and gga-miR-lc 
differ in their 3' ends and middle sequences by one nucleo- 
tide, while gga-miR-133a and gga-miR-133c differ in their 
3' end sequences by only one nucleotide. These sequence 
polymorphisms result in differential cleavage of the same 
miRNA precursor to produce a larger number of miRNA 
target sites. Gga-miR-la-1 and miR-la-2 share identical ma- 
ture sequences but are located on different chromosomes. 
Many mRNAs can be targets of one miRNA, and different 
combinations of miRNAs may coordinate to regulate spe- 
cific target genes [37]. For example, miR-l and miR-206 
have four mutual targets: histone deacetylase 4 (HDAC4) 
[38, 39], connexin43 (Cx43) [40], paired box 7 (Pax7 ) [41], 
and cMet [42]. Additionally, miR-206 can control myocyte 
differentiation by inhibiting fibroblast suppressors such as 
follistatin- 1 (Fstll), utrophin (Utrn) [43], and a subunit of 
DNA polymerase alpha (Polal) [44], indicating its extensive 
role the complicated gene network. 

Some miRNAs are expressed ubiquitously in different 
tissues, while some are expressed in a tissue-specific man- 
ner; miR-l, miR-133, and miR-206 belong to the latter cate- 
gory. These muscle-specific miRNAs have been demon- 
strated to regulate skeletal development in different manners 
[38, 45] and were dramatically up-regulated during myoblast 
differentiation [39, 46]. Up-regulation of miR-l is primarily 
due to accumulation of myocyte enhancer factor-2 (MEF2) 
resulting from decreased expression of HDAC4; interest- 
ingly, this confined expression of HDAC4 is caused by high 
miR-l expression. Thus, when MEF2 expression increases 
accompanied by up-regulation of miR-l, HDAC4 activity is 
repressed, resulting in further expression of MEF2 and con- 
trol of myocyte differentiation. 

3.2. Identification and Characterization of miRNAs Af- 
fecting Chicken Growth 

Previous studies have been focused on identifying and 
characterizing a large number of miRNAs associated with 
embryo development [47, 48], skeletal muscle development 
[49], lipogenesis and cell proliferation [50], lung and trachea 
with avian influenza virus infection [51], and embryo fibro- 
blasts infected with Marek's disease [52]. In the absence of 
differentially expressed miRNAs, gene expression would be 
dominated by monotony at the post-transcriptional and transla- 
tional levels, and the resulting message would be unreadable 
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and meaningless. Differential miRNAs found at different time 
points in the developing embryo and in discriminating skeletal 
muscle are of significant concern with regard to chicken 
growth. Generally, skeletal muscle goes through two stages of 
myoblast proliferation and differentiation during the embry- 
onic phase and postnatal muscle fiber hypertrophy. During the 
embryonic phase, the appearance of somites represents muscle 
development. Progenitor cells derived from somites generate 
myoblasts, which differentiate into mature muscle fibers via 
proliferation, migration, and fusion [53]. Muscle growth after 
birth is accomplished primarily through resizing of muscle 
fibers resulting from satellite cell fusion to existing fibers, 
rather than changing their numbers, as the final number of 
muscle fibers is determined at the end of embryogenesis [54]. 

As described above, miR-1 and miR-133 are involved in 
skeletal muscle growth and development. A bicistronic gene 
cluster encoding miR-1 and miR-133a, and another encoding 
miR-133b and miR-206 are transcribed from non-coding re- 
gions on mouse chromosomes 2 and 1, respectively. Tran- 
scription of these miRNAs is mainly achieved by myogenic 
transcription factors, such as myogenic differentiation 1 
(MyoD), serum response factor (SRF), and MEF2 by acting 
through cis-regulatory elements [55-57]. For instance, tran- 
scription regulation of miR-111331206 can be implemented by 
mammalian target of rapamycin (mTOR) signaling in a Myo- 
dependent manner. With the aid of MyoD, which lies down- 
stream of mTOR, transcription of miR-1 is triggered and in 
turn degrades HDAC4, a follistatin suppressor, affecting myo- 
cyte fusion. Similarly, inhibiting mTOR by rapamycin elicits a 
sharp decrease in miR-133 and miR-206 [58]. Moreover, miR- 
1 and miR-133 can regulate myogenesis via phosphatidylinosi- 
tol 3-kinase/v-AKT murine thymoma viral oncogene (Akt) 
(PDK/Akt) signaling Fig. (1), which at least partially partici- 
pates in muscle growth and hypertrophy [59]. Through IGF-1 
stimulation in mouse neonatal cardiomyocytes (CMCs) and 
during C2C12 skeletal muscle cell differentiation, binding of 
IGF-1 and IGF-1R results in activation of the PDK/Akt path- 
ways via PI3K-dependent phosphorylation and activation of 
Akt [60], stimulating hypertrophic growth and glucose uptake. 
In turn, active Akt phosphorylates and depresses transcription 
factor forkhead box 03 (Foxo3a), a negative regulator of pro- 
tein synthesis and muscle growth which regulates miR-1 pro- 
moter activity; the declining miR-1 reciprocally regulates IGF- 
1 and IGF-1R expression since they are targets of miR-1. This 
inverse correlation between IGF-1 protein levels and miR-1 is 
also found in myocardial biopsies of acromegalic patients 
[61]. Moreover, inverse regulation of miR-133 and IGF-1R 
was observed during C2C12 skeletal muscle cell differentia- 
tion [62]. Expression of miR-133 is up-regulated in response 
to increased expression of myogenin, a myogenic transcription 
factor which transactivates miR-133 in the presence of exoge- 
nous IGF-1. Overexpression of miR-133 has important impli- 
cations on the IGF-IR/PDK/Akt signaling pathway due to its 
ability to block IGF-1R expression and sequentially decrease 
regulation of Akt phosphorylation, miR-1 overexpression 
gives rise to pronounced down-regulation of phospho-Akt 
supported by a generalized decrease in the IGF- 1 signal trans- 
duction pathway. Another muscle-enriched miRNA, miR-486, 
reportedly indirectly and positively inhibits activation of 
PI3K/Akt signaling in rat CMCs through translation repres- 
sion of two crucial negative regulators, phosphatase and tensin 



homolog (PTEN) and Foxola [63]. If activated by MyoD and 
my ocardin -related transcription factor-A (MRTF-A), accumu- 
lation of miR-486 will rise and activate PI3K, target PTEN, or 
activate the phosphorylation of Akt by targeting Foxola to 
activate pathway. Active Akt in reverse inhibits Foxola activ- 
ity. Thus, miR-1 and miR-206 play important roles in promot- 
ing differentiation. In contrast, miR-133 overexpression inhibits 
myoblast differentiation, but promotes myoblast proliferation by 
targeting SRF, a regulator of miR-133 transcription [39]. 




Myogenin 



SRF. MRTF-A. MyoD 



MyoD, myogenin. MEF2, SRF 



Fig. (1). MicroRNAs regulating myogenesis via PDK/Akt signal- 
ing. This schematic diagram is a reference for [60-63]. The dash 
arrows indicates the favorable pathway has not been well clarified. 
IGF-1: insulin-like growth factor 1; IGF-1R: insulin-like growth 
factor 1 receptor; AKT P : phospho-AKT; Foxola: forkhead box 01; 
Foxo3a: forkhead box 03; MEF2: myocyte enhancer factor-2; 
MRTF-A: myocardin-related transcription factor-A; SRF: serum 
response factor; PTEN: phosphatase and tensin homolog; PI3K: 
phosphoinositide-3-kinase; MyoD: myogenic differentiation 1. 

3.3. MicroRNAs and Embryonic Development 

Dynamic sets of miRNAs have been demonstrated to 
express, which influence embryogenesis and organogenesis 
in chicks at 11 days of incubation [48]. By constructing and 
sequencing a small RNA library, a many known miRNAs 
(approximately 38% of the embryonic chick small RNA li- 
brary), along with small regions of homologous miRNAs 
with other species and novel miRNAs, were identified. Iden- 
tified miRNAs such as miR-143, miR-214, miR-363, miR- 
10a, and miR-22 show a significant degree of conservation 
compared to mouse, gorilla, bovine, and human. At least 27 
existing miRNA clusters were discovered in the chicken ge- 
nome, and nearly all are conserved within vertebrate species, 
suggesting that miRNAs both have an evolutionally ancient 
origin and share similar roles among species [64]. Novel 
miRNAs represented 1% of identified miRNAs, while miR- 
125b accounted for 11% of the total library [48]. The most 
abundant miR-125b is significant. Most recently, miR-125b 
was implicated in myoblast differentiation in vitro and mus- 
cle regeneration in vivo. Expression of IGF-2 can be directly 
regulated by enhanced mTOR signaling at the transcriptional 
level [65]. A recent study complemented miR-1 25b to corre- 
late IGF-2 and mTOR [66], i.e., expression of miR-125b is 
negatively controlled by mTOR signaling in an mTOR 
kinase-independent manner. When mTOR signaling is en- 
hanced, miR-1 25b is reduced, and the inhibitory effect of 
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miR-125b on IGF-2 is weakened, leading to high IGF-2 ex- 
pression. IGF-2 is an autocrine factor that initiates myoblast 
differentiation in vitro [67]; its high expression is imperative 
for injury-induced muscle regeneration [68]. 

In situ hybridization of chicken embryos during early 
stages of incubation revealed that miR-1, which facilitates 
differentiation of mesodermal progenitors to the muscle line- 
age [69], was detected in the somatic myotome beginning at 
stage 14 [70], which corresponds to the onset of skeletal 
muscle cell differentiation. In contrast, miR-133a was de- 
tected in the myocardium and myotome at stage 15, coinci- 
dent with rapid progress in limb development, and miR-206 
was detected in somites at stage 20, which is characterized 
by limb-bud enlargement [71]. Furthermore, restricted ex- 
pression of miR-206 steadily increased over 1.5 to 5 d of 
incubation and was exclusive in developing somites, particu- 
larly in the developing myotome. An opposite trend in ex- 
pression was observed following fibroblast growth factor 
(FGF) bead implantation due to FGF-mediated signaling, 
which negatively regulated the onset of miR-206 expression 
[72]. FGF signals can appoint some sets of progenitor cells 
in the ventral somite to ribs and tendons by acting through 
the mitogen-activated protein kinase (MAPK) signal trans- 
duction cascade [73]. Temporal expression of specific 
miRNA is necessary for organ development; therefore, 
growth and development of the embryo can effectively oc- 
cur. However, during embryogenesis and organogenesis, 
formation of the somite is a key step. 

Somites are short-lived mesodermal structures for which 
ventral somite cells engender the sclerotome containing pro- 
genitor cells for cartilage and bone, and dorsal somite cells 
form the dermomyotome containing progenitor cells for 
skeletal muscle. This acts as a visible symbol for the seg- 
mented nature of the vertebrate body plan [74]. Individual 
vertebrae are derived from the somite, which is generated 
following sequential segmentation of a region of the preso- 
matic mesoderm. Because few novel miRNAs were identi- 
fied in somites using traditional sequencing, Solexa sequenc- 
ing was used to globally explore novel miRNAs involving 
chicken somites [75]. Conventional sequencing of cDNA is a 
cumbersome process and cannot easily detect miRNAs pre- 
sent in relatively low abundance; deep sequencing or mas- 
sively parallel signature sequencing (MPSS) can be used to 
overcome these limitations [76]. Rathjen et al. obtained 
651,273 reads from dissected somite tissue from 3-, 4-, and 
5-day old embryos and identified 42 new miRNA candidates 
[75]. Eighteen were confirmed using Northern blot analysis. 
Notably, the novel miR-lOa was confirmed with a high num- 
ber of reads (28,660), and miR-lOb was the most abundant 
(113,106) among known miRNAs, implying that these two 
miRNAs are related during somite development. Effects 
caused by miRNAs on somite development and embryonic 
morphology have been reported previously, such as those 
caused by miR-196 [77]. miR-196 acts upstream of the ho- 
meobox (Hox) gene to define the boundary of Hox gene ex- 
pression along the anterior-posterior embryonic axis. A cen- 
tral function of the spatio-temporal expression of Hox gene 
in embryonic development is to pattern of the vertebrate ax- 
ial skeleton. Unilaterally localized injection of an- 
tagomiR196 into the presomatic mesoderm in ova leads to 
vertebral transformations at the cervical-thoracic boundary, 
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and thus high-frequency ectopia of the last cervical vertebra 
(cl4) occurs as decreased miR-196 expression gives expands 
the limit of homeoboxB8 (Hoxb8) transcripts in more ante- 
rior somites. Intriguingly, miR-10 has been shown decreases 
expression of several Hox genes in zebrafish [78]. Whether 
this data agrees with the mechanism of regulation of Hox 
targets with miR-196 requires further investigation since 
little evidence for its role in embryogenesis has been found. 

Additionally, Glazov et al. found that the total number of 
let-7 family reads increased significantly during embryonic 
stages on days 5, 7, and 9, particularly let-7b, with maximum 
reads on day 9 [47]. After identification in C. elegans, let-7 
miRNAs have been shown to play vital roles in mediating 
cell proliferation and differentiation. As a lung tumor sup- 
pressor in humans, overexpression of let-7 results in inverse 
expression of RAS, a potential oncogene, and inhibits lung 
tumor cell growth [79]. It was reported to exert a negative 
effect on cell number and positive effect on the fraction of 
cells in the G2/M cell cycle phase following peripheral in- 
troduction of let-7 in primary fibroblasts by targeting and 
down-regulating the cell division cycle 34 (Cdc34) gene, 
indicating its crucial influence on cell cycle control [80]. 
Introduction of let-7 family members in mouse embryonic 
stem cells (ESCs) can suppress continuous self-renewal re- 
sulting from a lack of DiGeorge syndrome critical region 
gene 8 (Dgcr8), which enables silencing of this program. 
Moreover, inhibition of let-7b boosts de-differentiation of 
somatic cells into induced pluripotent stem (iPS) cells [81]. 
Overexpression of let-7b elevates mouse neural stem cell 
differentiation but reduces proliferation by targeting the stem 
cell regulator nuclear receptor subfamily 2, group E, mem- 
ber 1 (Nr2el) and the cell cycle regulator cyclin D7[82]. 
This information underscores its specific role in stabilizing 
the self -renewing program and differentiating stem cells. 
Two miRNAs, miR-1623 and miR-181b, are differentially 
expressed in 14-day-old embryos of diverse growth-rate 
dwarf and normal chickens [83]. 

In summary, miRNAs, such as miR-101,p\ay an indispen- 
sable role in establishing certain organs by localizing or fast 
shifting expression in appropriate timing with embryonic de- 
velopment. miR-101 may monitor gonadal development to 
determine the gender of a chicken [84]. Expression increases 
between 5.5 d and 9.5 d of the embryonic phase in both sexes, 
with higher expression in males at E5.5 and E7.5, but a sig- 
nificant increase in females is observed at E9.5. It is likely that 
the high abundance of miR-101 after gonadal differentiation in 
females is crucial for determining the nature of ovarian cells 
due to its inhibitory effect on SRY ( sex determining region Y)- 
box 9 (SOX9), a key component of testes differentiation. Addi- 
tionally, upon suppression of some signaling inhibitory fac- 
tors, transforming growth factor beta/ Anti-Mullerian hormone 
(TGF-[3/AMH) signaling is reinforced to regulate follicle acti- 
vation. In males, miR-101 may be responsible for testes forma- 
tion by fine tuning SOX9. As a sensitive indicator, spatio- 
temporal expression of miRNAs reflects the accurate devel- 
opmental conditions of some organs. 

3.4. MicroRNAs and Growth of Skeletal Muscle 

Wang et al. identified 32 known miRNAs from skeletal 
muscle of Arbor Acres commercial chickens, of which 12 
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form five clusters: miR-133a-l-miR-la-2, miR-23b-miR-24, 
miR-99a-let-7c, miR-92-miR-19b-miR-18a-miR-17, and miR- 
30e-miR-30c-l, suggesting that most miRNAs co-express in 
skeletal muscle [85]. Because of the complexity of muscle 
development based upon spatiotemporal networks regulating 
mRNA transcription and translation, miRNA expression 
during eight different developmental stages (from the birth to 
7 w of age) was determined by using quantitative RT-PCR. 
Both miR-133a and miR-la were more important during later 
stages of development than during early stages of develop- 
ment since expression of both increased from 14 to 49 d and 
decreased from 0 to 14 d. 

Divergent muscle growth rate is the most common char- 
acteristic in commercial chicken lines, particularly in broilers 
and layers, normal, and SLD chickens. Broilers are raised for 
meat while layers are used for egg production, which is a 
consequence of the higher growth rate and larger muscle 
mass as well as bigger body size of broilers compared to 
layers. Both exterior selection pressure and intrinsic genetic 
patterns collectively elicit lager batteries of differences be- 
tween these two breeds. For SLD chickens, aberrant muscle 
development attributed to a defective GHR gene show typi- 
cal features of dwarf, serious weight loss, and decreased 
number of muscle fibers and fiber diameter. 

A recent preliminary study examined the miRNA tran- 
scriptome to identify differentially expressed miRNAs in- 
volved in skeletal muscle growth between broilers and layers 
[49]. A total of 33 novel chicken miRNAs were obtained, and 
102 miRNAs showed significant differences between broilers 
and layers. To confirm the accuracy of these results, 17 of 
these miRNAs were randomly validated using microarray 
analysis as well as real-time RT-PCR. Expression patterns of 
15 miRNAs identified using RT-PCR agreed with those iden- 
tified using deep sequencing, miR-101, miR-lOa, miR-lOb, 
miR-1677, let-7f, and miR-31 were higher in layers, while 
miR-200b, let-7c, miR-16c, miR15b, miR-15c, miR460, miR- 
429, miR-2188, and the novel miR-N2 were higher in broilers. 
Of note, miR-206 was the most abundant miRNA in both 
broilers and layers (131,609 reads and 222,998 reads, respec- 
tively), but significant differences in expression were ob- 
served. Expression of miR-133 was lower than that of miR-206 
and showed significant differences, miR-la expression was 
high in both libraries and showed significant differences in 
expression as well. Based on this information, the correspond- 
ing mRNA transcriptome was investigated and 57 candidate 
targets for 16 miRNAs were identified. Additionally, the final 
integrated network was constructed to depict miRNA-protein 
interactions and protein-protein interactions. Although three 
miRNAs, miR-1, miR-101, and miR-499, were predicted to 
target the activin A receptor type IIB (ACVR2B) gene, valida- 
tion was only performed for miR-1 and ACVR2B. We exam- 
ined miRNA expression in fast-growing White Recessive 
Rock chickens in contrast to slow-growing Qingyuan Par- 
tridge chickens at 7 w of age. We identified 23 differentially 
expressed miRNAs in the breast muscle of slow-growing 
chickens compared to in fast-growing chicken Fig. (2A). Wnt 
and insulin signaling pathways were two main pathways ac- 
cording to pathway analysis (data not published). 

The sex-linked dwarfism chicken is characterized by its 
relatively small body size, lower feed intake, lower basal me- 



tabolism, reduced IGF-1 levels, and high concentrations of 
GH in the plasma compared to in normal chickens. Previous 
studies suggested that mutations in the GHR gene may result 
in phenotypic variation of SLD [86-88]. These mutations in- 
duced either different transcripts in terms of transcription level 
or chaotic GHR proteins in terms of translation level. Com- 
parison of mRNA expression for GH, GHR, and IGF-1 genes 
at 56 d of age in the livers of sex-linked dwarf chickens with 
normal chickens showed that GH expression between these 
chickens showed equally high transcription levels. However, 
expression of GHR in dwarf chickens was significantly higher 
than that in normal chickens, and very little IGF-1 gene ex- 
pression was observed in dwarf chickens, further illustrating 
that the GH-independent dwarfism is due to dysfunctional 
GHR protein in accordance with a mutant GHR gene [89]. 
Notably, GHR-deficient mice showed reductions in myofiber 
number, coupling with defective muscle skeletal development 
due to diminished myoblast fusion [90]. Based upon our pre- 
vious study, miRNA regulation, which implicates GHR as 
well, may be important in dwarfism formation [83]. A signifi- 
cant difference in expression of let-7b in skeletal muscle was 
observed between SLD chickens due to the identical mutant 
type of GHR with Connecticut SLD [88] and normal chickens 
at 7 w of age. Comparison of differential expression profiles 
of miRNA and mRNA in SLD and normal chickens showed 
that, in addition to let-7b, mRNA expression of GHR genes 
increased SLD compared to in normal chicken. Because aber- 
rant GHR elicits the dwarfism phenotype, it is important in the 
functional mechanism that how GHR acts on skeletal muscle 
in SLD in response to let-7b. The relationship between let-7b 
and GHR has been confirmed through a dual reporter assay 
and overexpression experiment in vitro. As expected, the 
complementary sequence of let-7b to normal GHR matched 
the deleted region of mutant GHR, indicating that let-7b sup- 
presses expression of GHR in normal chickens but not in SLD 
chickens. The primary pathway in which GHR participates is 
for the janus kinase/signal transducers and activators of tran- 
scription (JAK/STAT) signaling pathway. Interestingly, 
mRNA expression of suppressor of cytokine signaling 3 
(SOCS3) was markedly up-regulated in SLD in contrast to 
other genes involved in this pathway. Thus, tight control of 
let-7b may block expression of SOCS3 in a non-mutant GHR- 
dependent manner. 

SOCS3 is a suppressor of cytokine signaling which blocks 
the JAK/STAT pathway relying on binding to janus kinase 
and to specific cytokine receptors, which affects cell prolifera- 
tion, differentiation, apoptosis, and immunoregulation [91]. 
SOCS3 also blocks insulin signaling by targeting insulin re- 
ceptor substrate 1 (IRS1) and insulin receptor substrate 2 
(IRS2), two critical signaling molecules for insulin action, and 
forces their ubiquitination and degradation, contributing to 
insulin resistance [92] . Leptin typically influences the feeding 
and neuroendocrine function depending on binding to and 
consequently activating its receptor (LEPR). For SLD, the 
indirect, acute up-regulation of SOCS3 impairs the effect of 
leptin by inhibiting LEPR production. Thus, leptin suppression 
increases feeding and glucocorticoid production, but decreases 
energy expenditure of reproductive, growth, and metabolic 
rate. SLD chickens eat less than normal chickens, perhaps for 
because of their smaller somatotypes or due to unknown regu- 
lators compensating for the effect induced by SOCS3. How- 
ever, the data provide insight into the regulatory mechanism. 
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Fig. (2). Differentially expressed miRNAs in chickens with different growth performance and target validation. A: Expression pattern of 23 
miRNAs in breast muscles of White Recessive Rock (WRR) and Qingyuan Partridge chickens (QYP) based on microarray experiments. The 
bars represent the fold change in WRR/QYP. B: Relative abundance of 22 miRNAs in the livers of normal and RSS chickens based on Solexa 
sequencing data. The bars represent the fold change of RSS/normal. C: Schema of gga-miR-221 binding sites in chicken wide-type and mu- 
tated TAB2 (TGF-beta activated kinase 1/MAP3K7 binding protein 2) 3'-UTR sequence. Seed binding sites and mutated bases are high- 
lighted in red. Target prediction was performed using TargetScan 5.0 software (http://targetscan.org/). D: To verify whether gga-miR-221 
targets TAB2, a luciferase reporter assay was conducted. A gga-miR-221 overexpression plasmid was constructed using pcDNA3.1 (Promega, 
Fitchburg, WI, USA). Wild-type (WT) or mutated TAB2 3'-UTR sequence (MUT) was cloned into the pmirGLO luciferase reporter vector at 
downstream of Renilla luciferase, the internal Firefly luciferase gene was used to normalize for transfection efficiency. DF-1 cells were used 
for transfection, and 600 ng of the gga-miR-221 overexpression plasmid was co-transfected with 200 ng of WT or MUT reporter plasmid 
(221-WT/-MUT). A parallel inhibition or negative control experiment was conducted using30 ng IN or NC vector plasmid cotransfected with 
200 ng of WT reporter plasmid (221-IN/-NC). Dual-luciferase assays were conducted in triplicate in 24-well plates. Cells were lysed and 
assayed for luciferase activity 36h after transfection. Error bars represent the standard error of three independent experiments per group. 



Notably, disease in chicken populations can also give rise 
to a dramatic drop in skeletal muscle mass and result in runt- 
ing and stunting syndrome (RSS). This is one of the most 
common syndromes causing serious growth obstruction, and 
chickens with this syndrome consume large amounts of feed, 
show poor growth performance, and can be easily affected 
by a disease. Thus far, its aetiology remains unclear. To 
identify the extensive role of miRNAs in balancing this 
complicated pathological mechanism, Solexa sequencing 
was performed to determine differential expression profiles 
of 22 miRNAs in the livers of normal and RSS chickens at 7 
w of age Fig. (2B). Subsequently, we found that TGF-beta 



activated kinase 1/MAP3K7 binding protein 2 (TAB2) was 
target of gga-miR-221 using a luciferase reporter gene assay 
Fig. (2C and 2D). Interestingly, miR-221 affected skeletal 
muscle cell growth and disease. It was down-regulated dur- 
ing myogenic differentiation, with assistance from the RAS- 
mitogen-activated protein kinase (MAPK) pathway. Ectopic 
expression of miR-221 in differentiated myoblasts extended 
the cell cycle and delayed myogenin expression, which had a 
suppressive role on cell cycle inhibitor p27 [93]. In chronic 
myeloid leukemia (CML) patients, miR-221was significantly 
up-regulated during the blast crisis phase and was predicted 
to target v-crk sarcoma virus CT10 oncogene homolog 
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(avian)-like (CRKL) in the MAPK signaling and phosphoi- 
nositide-3-kinase, regulatory subunit 1 (PIK3R1) in epider- 
mal growth factor receptor (EGFR) signaling associated with 
CML [94] . The roles of these two miRNAs and their targets 
underlying RSS must be further examined. In general, miR- 
NAs are thought to be involved in regulating multiple genes 
regardless of consistently expressing or differentially ex- 
pressing in different skeletal muscle by targeting and inhibit- 
ing their mRNAs. Thus, it is necessary to examine the func- 
tion and characteristics of each miRNA. 

4. CONCLUSION AND PROSPECTS 

Various studies have been conducted to identify connec- 
tions between genotype and phenotype. Chicken growth is 
complex and orchestrated through precise control of multiple 
genes. Association analysis using candidate gene and ge- 
nome scan using linkage analysis are two primary means of 
detecting QTL related to growth, but this method exhibits 
some limitations. Because candidate genes are limited when 
using candidate gene approaches and QTL mapping is im- 
precise due to insufficient marker density when using link- 
age analysis, major genes leading to phenotype differences 
can not be determined. With the development of biotechnol- 
ogy and bioinformatic approaches, GWAS can be used to 
identify genes using tens of thousands of SNP markers at the 
whole-genome level. GGA1 and GGA4 are two critical 
chromosomes affecting chicken growth, particularly body 
weight. As integral members of gene networks, miRNAs 
explain the complexity of life. Many growth-related miR- 
NAs have been discovered, including miR-1, miR-133, miR- 
206, miR-101, and let-7b, the biochemical roles of which 
have been demonstrated through experimental validation. In- 
depth knowledge of miRNAs will increase the understanding 
of the animals' growth molecular mechanisms. 

CONFLICT OF INTEREST 

The authors confirm that this article content has no con- 
flicts of interest. 

ACKNOWLEDGEMENTS 

This work was supported by the China Agriculture Re- 
search System (CARS-42-G05), China High-Tech Programs 
(2011AA100301), and Natural Scientific Foundation of 
China (31172200). 

REFERENCES 

[11 Natt, D.; Rubin, C.J.; Wright, D.; Johnsson, M; Belteky, J.; 
Andersson, L.; Jensen, P. Heritable genome-wide variation of gene 
expression and promoter methylation between wild and domesti- 
cated chickens. BMC Genomics, 2012, 13, 59. 

[2] Franklin, T.B.; Mansuy, I.M. Epigenetic inheritance in mammals: 
evidence for the impact of adverse environmental effects. Neuro- 
biolDis., 2010, 39, 61-65. 

[3] Molinier, J.; Ries, G.; Zipfel, C; Hohn, B. Transgeneration mem- 
ory of stress in plants. Nature, 2006, 442, 1046-1049. 

[41 Sato, S.; Ohtake, T.; Uemoto, Y.; Okumura, Y.; Kobayashi, E. 
Polymorphism of insulin-like growth factor 1 gene is associated 
with breast muscle yields in chickens. Anim Sci J, 2012, 83, 1-6. 

[5] Song, C; Gu, X.; Feng, C; Wang, Y.; Gao, Y.; Hu, X.; Li, N. 
Evaluation of SNPs in the chicken HMGA2 gene as markers for 
body weight gain. Anim. Genet., 2011, 42, 333-336. 

[6] Nie, Q.H.; Fang, M.X.; Xie, L.; Shen, X.; Liu, J.; Luo, Z.P.; Shi, 
J.J.; Zhang, X.Q. Associations of ATGL gene polymorphisms with 
chicken growth and fat traits. JAppl Genet, 2010, 51, 185-191. 



[7] Fang, M.; Nie, Q.; Luo, C; Zhang, D.; Zhang, X. Associations of 
GHSR gene polymorphisms with chicken growth and carcass traits. 
Mot Biol. Rep., 2010, 37, 423-428. 

[8] Nie, Q.; Fang, M.; Xie, L.; Zhou, M.; Liang, Z.; Luo, Z.; Wang, G.; 
Bi, W.; Liang, C; Zhang, W.; Zhang, X. The PIT1 gene polymor- 
phisms were associated with chicken growth traits. BMC Genet, 
2008, 9, 20. 

[9] Lei, M.; Peng, X.; Zhou, M.; Luo, C; Nie, Q.; Zhang, X. Polymor- 
phisms of the IGF1R gene and their genetic effects on chicken 
early growth and carcass traits. BMC Genet, 2008, 9, 70. 

[10] Ouyang, J.H.; Xie, L.; Nie, Q.; Luo, C; Liang, Y.; Zeng, H.; 
Zhang, X. Single nucleotide polymorphism (SNP) at the GHR gene 
and its associations with chicken growth and fat deposition traits. 
Br Poult Sci, 2008, 49, 87-95. 

[11] Li, C.C.; Li, K; Li, J; Mo, D.L.; Xu, R.F.; Chen, G.H.; Qiangba, 
Y.Z.; Ji, S.L.; Tang, X.H.; Fan, B.; Zhu, M.J.; Xiong, T.A.; Guan, 
X.; Liu, B. Polymorphism of ghrelin gene in twelve Chinese in- 
digenous chicken breeds and its relationship with chicken growth 
traits. Asian-Aust. J. Anim, 2006, 153-159. 

[12] Lei, M; Luo, C; Peng, X.; Fang, M.; Nie, Q.; Zhang, D.; Yang, G.; 
Zhang, X. Polymorphism of growth-correlated genes associated 
with fatness and muscle fiber traits in chickens. Poult Sci, 2007, 86, 
835-842. 

[13] Fang, M.; Nie, Q.; Luo, C; Zhang, D.; Zhang, X. An 8bp indel in 
exon 1 of Ghrelin gene associated with chicken growth. Domest 
Anim Endocrinol, 2007, 32, 216-225. 

[14] Qiu, F.F.; Nie, Q.H.; Luo, C.L.; Zhang, D.X.; Lin, S.M.; Zhang, 
X.Q. Association of single nucleotide polymorphisms of the insulin 
gene with chicken early growth and fat deposition. Poult Sci, 2006, 
55,980-985. 

[15] Nie, Q.; Sun, B.; Zhang, D.; Luo, C; Ishag, N.A.; Lei, M.; Yang, 
G.; Zhang, X. High diversity of the chicken growth hormone gene 
and effects on growth and carcass traits. /. Hered., 2005, 96, 698- 
703. 

[16] Lei, M.M.; Nie, Q.H.; Peng, X.; Zhang, D.X.; Zhang, X.Q. Single 
nucleotide polymorphisms of the chicken insulin-like factor bind- 
ing protein 2 gene associated with chicken growth and carcass 
traits. Poult Sci, 2005, 84, 1191-1198. 

[17] Amills, M.; Jimenez, N.; Villalba, D.; Tor, M.; Molina, E.; Cubilo, 
D.; Marcos, C; Francesch, A.; Sanchez, A.; Estany, J. Identifica- 
tion of three single nucleotide polymorphisms in the chicken insu- 
lin-like growth factor 1 and 2 genes and their associations with 
growth and feeding traits. Poult Sci, 2003, 82, 1485-1493. 

[18] Ankra-Badu, G.A.; Shriner, D.; Le Bihan-Duval, E.; Mignon- 
Grasteau, S.; Pitel, F.; Beaumont, C; Duclos, M.J.; Simon, J.; Por- 
ter, T.E.; Vignal, A.; Cogburn, LA.; Allison, D.B.; Yi, N; Aggrey, 
S.E. Mapping main, epistatic and sex-specific QTL for body com- 
position in a chicken population divergently selected for low or 
high growth rate. BMC Genomics, 2010, 11, 107. 

[19] Rao, Y.; Shen, X.; Xia, M.; Luo, C; Nie, Q.; Zhang, D.; Zhang, X. 
SNP mapping of QTL affecting growth and fatness on chicken 
GGA1. Genet. Set Evol, 2007, 39, 569-582. 

[20] Jennen, D.G.; Vereijken, A.L.; Bovenhuis, H.; Crooijmans, R.M.; 
van derPoel, J.J.; Groenen, MA. Confirmation of quantitative trait 
loci affecting fatness in chickens. Genet. Sel. Evol, 2005, 37, 215- 
228. 

[21] Uemoto, Y.; Sato, S.; Odawara, S.; Nokata, H.; Oyamada, Y.; 
Taguchi, Y.; Yanai, S.; Sasaki, O.; Takahashi, H.; Nirasawa, K.; 
Kobayashi, E. Genetic mapping of quantitative trait loci affecting 
growth and carcass traits in F2 intercross chickens. Poult Sci, 2009, 
88, 477-482. 

[22] Liu, X.; Li, H.; Wang, S.; Hu, X.; Gao, Y.; Wang, Q.; Li, N; 
Wang, Y.; Zhang, H. Mapping quantitative trait loci affecting body 
weight and abdominal fat weight on chicken chromosome one. 
Poult Sci, 2007, 86, 1084-1089. 

[23] Liu, X.; Zhang, H.; Li, H.; Li, N; Zhang, Y.; Zhang, Q.; Wang, S.; 
Wang, Q.; Wang, H. Fine-mapping quantitative trait loci for body 
weight and abdominal fat traits: effects of marker density and sam- 
ple size. Poult Sci, 2008, 87, 1314-1319. 

[24] Zhang, H.; Liu, S.H.; Zhang, Q.; Zhang, Y.D.; Wang, S.Z.; Wang, 
Q.G.; Wang, Y.X.; Tang, Z.Q.; Li, H. Fine-mapping of quantitative 
trait loci for body weight and bone traits and positional cloning of 
the RBI gene in chicken. J. Anim. Breed. Genet, 2011, 128, 366- 
375. 

[25] Xie, L.; Luo, C; Zhang, C; Zhang, R.; Tang, J.; Nie, Q.; Ma, L.; 
Hu, X.; Li, N.; Da, Y.; Zhang, X. Genome-wide association study 



Chicken Growth Traits Determined Based on GWAS and microRNA Regulation 



Current Genomics, 2013, Vol. 14, No. 2 145 



identified a narrow chromosome 1 region associated with chicken 
growth traits. PLoS One, 2012, 7, e30910. 
[26] Carlborg, O.; Kerje, S.; Schutz, K.; Jacobsson, L.; Jensen, P.; [48] 
Andersson, L. A global search reveals epistatic interaction between 
QTL for early growth in the chicken. Genome Res., 2003, 13, 413- [49] 
421. 

[27] Gu, X.; Feng, C; Ma, L.; Song, C; Wang, Y.; Da, Y.; Li, H.; 
Chen, K.; Ye, S.; Ge, C; Hu, X.; Li, N. Genome-wide association 
study of body weight in chicken F2 resource population. PLoS One, [50] 
2011,6, e21872. 

[28] Johansson, A.M.; Pettersson, M.E.; Siegel, P.B.; Carlborg, O. Ge- 
nome-wide effects of long-term divergent selection. PLoS Genet, [51] 
2010, 6, el001188. 

[29] Sladek, R.; Rocheleau, G.; Rung, J.; Dina, C; Shen, L.; Serre, D.; 
Boutin, P.; Vincent, D.; Belisle, A.; Hadjadj, S.; Balkau, B.; Heude, 
B.; Charpentier, G.; Hudson, T.J.; Montpetit, A.; Pshezhetsky, 

A. V.; Prentki, M.; Posner, B.I.; Balding, D.J.; Meyre, D.; Poly- [52] 
chronakos, C; Froguel, P. A genome-wide association study identi- 
fies novel risk loci for type 2 diabetes. Nature, 2007, 445, 881-885. 

[30] Zeggini, E.; McCarthy, M.I. Identifying susceptibility variants for 

type 2 diabetes. Methods Mol Biol, 2007, 376, 235-250. [53] 

[31] Kamei, Y.; Miura, S.; Suzuki, M.; Kai, Y.; Mizukami, J.; Tanigu- 
chi, T.; Mochida, K.; Hata, T.; Matsuda, J.; Aburatani, H.; Nishino, 
I.; Ezaki, O. Skeletal muscle FOXOl (FKHR) transgenic mice have 
less skeletal muscle mass, down-regulated Type I (slow twitch/red [54] 
muscle) fiber genes, and impaired glycemic control. J. Biol. Chem., 
2004,279,41114-41123. [55] 

[32] Filleur, S.; Hirsch, J.; Wille, A.; Schon, M.; Sell, C; Shearer, M.H.; 
Nelius, T.; Wieland, I. INTS6/DICE1 inhibits growth of human an- 
drogen-independent prostate cancer cells by altering the cell cycle 
profile and Wnt signaling. Cancer Cell Int, 2009, 9, 28. 

[33] Javerzat, S.; Franco, M.; Herbert, J.; Platonova, N.; Peille, A.L.; [56] 
Pantesco, V.; De Vos, J.; Assou, S.; Bicknell, R.; Bikfalvi, A.; 
Hagedorn, M. Correlating global gene regulation to angiogenesis in 
the developing chick extra-embryonic vascular system. PLoS One, [57] 
2009, 4, e7856. 

[34] le Sage, C; Agami, R. Immense promises for tiny molecules: un- 
covering miRNA functions. Cell Cycle, 2006, 5, 1415-1421. [58] 

[35] Filipowicz, W. RNAi: the nuts and bolts of the RISC machine. 
Cell, 2005, 122, 17-20. 

[36] Kim, V.N.; Han, J.; Siomi, M.C. Biogenesis of small RNAs in [59] 
animals. Nat Rev Mol Cell Biol, 2009, 10, 126-139. 

[37] Krek, A.; Grun, D.; Poy, M.N.; Wolf, R.; Rosenberg, L.; Epstein, [60] 
E.J.; MacMenamin, P.; Da, P.I.; Gunsalus, K.C.; Stoffel, M.; Ra- 
jewsky, N. Combinatorial microRNA target predictions. Nat. 
Genet., 2005, j7, 495-500. 

[38] Williams, A.H.; Liu, N.; van Rooij, E.; Olson, E.N. MicroRNA [61] 
control of muscle development and disease. Curr. Opin. Cell Biol, 
2009, 21, 461-469. 

[39] Chen, J.F.; Mandel, E.M.; Thomson, J.M.; Wu, Q.; Callis, T.E.; 
Hammond, S.M.; Conlon, F.L.; Wang, D.Z. The role of mi- 
croRNA- 1 and microRNA- 133 in skeletal muscle proliferation and 
differentiation. Nat. Genet, 2006, 38, 228-233. [62] 

[40] Anderson, C; Catoe, H.; Werner, R. MIR-206 regulates con- 
nexin43 expression during skeletal muscle development. Nucleic 
Acids Res., 2006, 34, 5863-5871. [63] 

[41] Chen, J.F.; Tao, Y.; Li, J.; Deng, Z.; Yan, Z.; Xiao, X.; Wang, D.Z. 
microRNA- 1 and microRNA-206 regulate skeletal muscle satellite 
cell proliferation and differentiation by repressing Pax7. /. Cell 
Biol., 2010, 790,867-879. [64] 

[42] Yan, D.; Dong, X.E.; Chen, X.; Wang, L.; Lu, C; Wang, J.; Qu, J.; 
Tu, L. MicroRNA- 1/206 targets c-Met and inhibits rhabdomyosar- 
coma development. J. Biol. Chem., 2009, 284, 29596-29604. [65] 

[43] Rosenberg, M.I.; Georges, S.A.; Asawachaicharn, A.; Analau, E.; 
Tapscott, S.J. MyoD inhibits Fstll and Utrn expression by inducing 
transcription of miR-206. J. Cell Biol., 2006, 175, 77-85. [66] 

[44] Kim, H.K.; Lee, Y.S.; Sivaprasad, U.; Malhotra, A.; Dutta, A. 

Muscle-specific microRNA miR-206 promotes muscle differentia- [67] 
tion. J. Cell Biol., 2006, 174, 677-687. 

[45] Chen, J.F.; Callis, T.E.; Wang, D.Z. microRNAs and muscle disor- 
ders. /. Cell Set, 2009, 122, 13-20. 

[46] Wong, C.F.; Tellam, R.L. MicroRNA-26a targets the histone meth- [68] 
yltransferase Enhancer of Zeste homolog 2 during myogenesis. J. 
Biol. Chem., 2008, 25^,9836-9843. 

[47] Glazov, E.A.; Cottee, PA.; Barris, W.C.; Moore, R.J.; Dalrymple, 

B. P.; Tizard, M.L. A microRNA catalog of the developing chicken 



embryo identified by a deep sequencing approach. Genome Res., 
2008,78,957-964. 

Hicks, J.A.; Tembhurne, P.; Liu, H.C. MicroRNA expression in 
chicken embryos. Poult Sci, 2008, 87, 2335-2343. 
Li, T.; Wu, R.; Zhang, Y.; Zhu, D. A systematic analysis of the 
skeletal muscle miRNA transcriptome of chicken varieties with di- 
vergent skeletal muscle growth identifies novel miRNAs and dif- 
ferentially expressed miRNAs. BMC Genomics, 2011, 72, 186. 
Hicks, J. A.; Trakooljul, N.; Liu, H.C. Discovery of chicken mi- 
croRNAs associated with lipogenesis and cell proliferation. 
Physiol. Genomics, 2010, 41, 185-193. 

Wang, Y.; Brahmakshatriya, V.; Zhu, H.; Lupiani, B.; Reddy, 
S.M.; Yoon, B.J.; Gunaratne, P.H.; Kim, J.H.; Chen, R.; Wang, J.; 
Zhou, H. Identification of differentially expressed miRNAs in 
chicken lung and trachea with avian influenza virus infection by a 
deep sequencing approach. BMC Genomics, 2009, 10, 512. 
Burnside, J.; Ouyang, M.; Anderson, A.; Bernberg, E.; Lu, C; 
Meyers, B.C.; Green, P.J.; Markis, M.; Isaacs, G.; Huang, E.; 
Morgan, R.W. Deep sequencing of chicken microRNAs. BMC Ge- 
nomics, 2008, 9, 185. 

Buckingham, M.; Bajard, L.; Chang, T.; Daubas, P.; Hadchouel, J.; 
Meilhac, S.; Montarras, D.; Rocancourt, D.; Relaix, F. The forma- 
tion of skeletal muscle: from somite to limb. J. Anat., 2003, 202, 
59-68. 

Moss, F.P.; Leblond, CP. Satellite cells as the source of nuclei in 
muscles of growing rats. Anat Rec, 1971, 770, 421-435. 
Liu, N.; Williams, A.H.; Kim, Y.; McAnally, J.; Bezprozvannaya, 
S.; Sutherland, L.B.; Richardson, J. A.; Bassel-Duby, R.; Olson, 
E.N. An intragenic MEF2-dependent enhancer directs muscle- 
specific expression of microRNAs 1 and 133. Proc Natl Acad Sci U 
SA, 2007, 704, 20844-20849. 

Rao, P.K.; Kumar, R.M.; Farkhondeh, M.; Baskerville, S.; Lodish, 
H.F. Myogenic factors that regulate expression of muscle-specific 
microRNAs. Proc Natl Acad Sci USA, 2006, 103, 8721-8726. 
Zhao, Y.; Samal, E.; Srivastava, D. Serum response factor regulates 
a muscle-specific microRNA that targets Hand2 during cardiogene- 
sis. Nature, 2005, 436, 214-220. 

Sun, Y.; Ge, Y.; Drnevich, J.; Zhao, Y.; Band, M.; Chen, J. Mam- 
malian target of rapamycin regulates miRNA- 1 and follistatin in 
skeletal myogenesis. J. Cell Biol, 2010, 789, 1 157-1 169. 
Bassel-Duby, R.; Olson, E.N. Signaling pathways in skeletal mus- 
cle remodeling. Annu. Rev. Biochem., 2006, 75, 19-37. 
Brunet, A.; Bonni, A.; Zigmond, M.J.; Lin, M.Z.; Juo, P.; Hu, L.S.; 
Anderson, M.J.; Arden, K.C.; Blenis, J.; Greenberg, M.E. Akt pro- 
motes cell survival by phosphorylating and inhibiting a Forkhead 
transcription factor. Cell, 1999, 96, 857-868. 

Elia, L.; Contu, R.; Quintavalle, M.; Varrone, F.; Chimenti, C; 
Russo, M.A.; Cimino, V.; De Marinis, L.; Frustaci, A.; Catalucci, 
D.; Condorelli, G. Reciprocal regulation of microRNA- 1 and insu- 
lin-like growth factor- 1 signal transduction cascade in cardiac and 
skeletal muscle in physiological and pathological conditions. Cir- 
culation, 2009, 720, 2377-2385. 

Huang, M.B.; Xu, H.; Xie, S.J.; Zhou, H.; Qu, L.H. Insulin-like 
growth factor-1 receptor is regulated by microRNA-133 during 
skeletal myogenesis. PLoS One, 2011, 6, e29173. 
Small, E.M.; O'Rourke, J.R.; Moresi, V.; Sutherland, L.B.; McA- 
nally, J.; Gerard, R.D.; Richardson, J. A.; Olson, E.N. Regulation of 
PI3-kinase/Akt signaling by muscle-enriched microRNA-486. Proc 
Natl Acad Sci V S A, 2010, 707,4218-4223. 

Shao, P.; Zhou, H.; Xiao, Z.D.; He, J.H.; Huang, M.B.; Chen, Y.Q.; 
Qu, L.H. Identification of novel chicken microRNAs and analysis 
of their genomic organization. Gene, 2008, 418, 34-40. 
Erbay, E.; Park, I.H.; Nuzzi, P.D.; Schoenherr, C.J.; Chen, J. IGF-II 
transcription in skeletal myogenesis is controlled by mTOR and nu- 
trients. J. Cell Biol, 2003, 163, 931-936. 

Ge, Y.; Sun, Y.; Chen, J. IGF-II is regulated by microRNA-125b in 
skeletal myogenesis. J. Cell Biol, 2011, 792, 69-81. 
Florini, J.R.; Magri, K.A.; Ewton, D.Z.; James, P.L.; Grindstaff, 
K; Rotwein, P.S. "Spontaneous" differentiation of skeletal 
myoblasts is dependent upon autocrine secretion of insulin-like 
growth factor-II. /. Biol. Chem., 1991, 266, 15917-15923. 
Paoni, N.F.; Peale, F.; Wang, F.; Errett-Baroncini, C; Steinmetz, 
H.; Toy, K; Bai, W.; Williams, P.M.; Bunting, S.; Gerritsen, M.E.; 
Powell-Braxton, L. Time course of skeletal muscle repair and gene 
expression following acute hind limb ischemia in mice. Physiol. 
Genomics, 2002, 77, 263-272. 



146 Current Genomics, 2013, Vol. 14, No. 2 



Xu et al. 



[69] Ivey, K.N.; Muth, A.; Arnold, J.; King, F.W.; Yeh, R.F.; Fish, J.E.; 
Hsiao, E.C.; Schwartz, R.J.; Conklin, B.R.; Bernstein, H.S.; Srivas- 
tava, D. MicroRNA regulation of cell lineages in mouse and human 
embryonic stem cells. Cell Stem Cell, 2008, 2, 219-229. 

[70] Hamburger, V.; Hamilton, H.L. A series of normal stages in the 
development of the chick embryo. 1951. Dev Dyn, 1992, 195, 231- 
272. 

[71] Darnell, D.K.; Kaur, S.; Stanislaw, S.; Konieczka, J.H.; 
Yatskievych, T.A.; Antin, P.B. MicroRNA expression during chick 
embryo development. Dev Dyn, 2006, 235, 3156-3165. 

[72] Sweetman, D.; Rathjen, T.; Jefferson, M.; Wheeler, G.; Smith, 
T.G.; Wheeler, G.N.; Munsterberg, A.; Dalmay, T. FGF-4 signal- 
ing is involved in mir-206 expression in developing somites of 
chicken embryos. Dev Dyn, 2006, 235, 2185-2191. 

[73] Smith, T.G.; Sweetman, D.; Patterson, M.; Keyse, S.M.; Munster- 
berg, A. Feedback interactions between MKP3 and ERK MAP 
kinase control scleraxis expression and the specification of rib pro- 
genitors in the developing chick somite. Development, 2005, 132, 
1305-1314. 

[74] Aulehla, A.; Herrmann, B.G. Segmentation in vertebrates: clock 
and gradient finally joined. Genes Dev, 2004, 18, 2060-2067. 

[75] Rathjen, T.; Pais, H.; Sweetman, D.; Moulton, V.; Munsterberg, A.; 
Dalmay, T. High throughput sequencing of microRNAs in chicken 
somites. Febs Lett., 2009, 583, 1422-1426. 

[76] Moxon, S.; Jing, R.; Szittya, G.; Schwach, F.; Rusholme, P.R.; 
Moulton, V.; Dalmay, T. Deep sequencing of tomato short RNAs 
identifies microRNAs targeting genes involved in fruit ripening. 
Genome Res., 2008, 18, 1602-1609. 

[77] McGlinn, E.; Yekta, S.; Mansfield, J.H.; Soutschek, J.; Bartel, D.P.; 
Tabin, C.J. In ovo application of antagomiRs indicates a role for miR- 
196 in patterning the chick axial skeleton through Hox gene regula- 
tion. Proc Natl Acad Sci USA, 2009, 106, 18610-18615. 

[78] Woltering, J.M.; Durston, A.J. MiR-10 represses HoxBla and 
HoxB3a in zebrafish. PLoS One, 2008, 3, el396. 

[79] Johnson, CD.; Esquela-Kerscher, A.; Stefani, G.; Byrom, M.; 
Kelnar, K.; Ovcharenko, D.; Wilson, M.; Wang, X.; Shelton, J.; 
Shingara, J.; Chin, L.; Brown, D.; Slack, F.J. The let-7 microRNA 
represses cell proliferation pathways in human cells. Cancer Res., 
2007,67,7713-7722. 

[80] Legesse-Miller, A.; Elemento, O.; Pfau, S.J.; Forman, J.J.; Tava- 
zoie, S.; Coller, H.A. let-7 Overexpression leads to an increased 
fraction of cells in G2/M, direct down-regulation of Cdc34, and 
stabilization of Weel kinase in primary fibroblasts. /. Biol. Chem., 
2009, 284, 6605-6609. 

[81] Melton, C; Judson, R.L.; Blelloch, R. Opposing microRNA fami- 
lies regulate self-renewal in mouse embryonic stem cells. Nature, 
2010,465,621-626. 

[82] Zhao, C; Sun, G.; Li, S.; Lang, M.F.; Yang, S.; Li, W.; Shi, Y. 



MicroRNA let-7b regulates neural stem cell proliferation and dif- 
ferentiation by targeting nuclear receptor TLX signaling. Proc Natl 
AcadSciUSA, 2010, 107, 1876-1881. 

[83] Lin, S.; Li, H.; Mu, H.; Luo, W.; Li, Y.; Jia, X.; Wang, S.; Jia, X.; 
Nie, Q.; Li, Y.; Zhang, X. Let-7b regulates the expression of the 
growth hormone receptor gene in deletion-type dwarf chickens. 
BMC Genomics, 2012, 13, 306. 

[84] Cutting, A.D.; Bannister, S.C.; Doran, T.J.; Sinclair, A.H.; Tizard, 
M.V.; Smith, C.A. The potential role of microRNAs in regulating 
gonadal sex differentiation in the chicken embryo. Chromosome 
Res., 2012, 20, 201-213. 

[85] Wang, X.G.; Yu, J.F.; Zhang, Y.; Gong, D.Q.; Gu, Z.L. Identifica- 
tion and characterization of microRNA from chicken adipose tissue 
and skeletal muscle. Poult Sci, 2012, 91, 139-149. 

[86] Hull, K.L.; Marsh, J. A.; Harvey, S. A missense mutation in the 
GHR gene of Cornell sex-linked dwarf chickens does not abolish 
serum GH binding. J. Endocrinol., 1999, 161, 495-501. 

[87] Tanaka, M.; Hayashida, Y.; Wakita, M.; Hoshino, S.; Nakashima, 
K. Expression of aberrantly spliced growth hormone receptor 
mRNA in the sex-linked dwarf chicken, Gifu 20. Growth Regul, 
1995,5,218-223. 

[88] Agarwal, S.K.; Cogburn, L.A.; Burnside, J. Dysfunctional growth 
hormone receptor in a strain of sex-linked dwarf chicken: evidence 
for a mutation in the intracellular domain. J. Endocrinol., 1994, 
142, 427-434. 

[89] Wu, G.Q.; Zheng, J.X.; Yang, N. [Expression profiling of GH, 
GHR, and IGF-1 genes in sex-linked dwarf chickens]. Yi Chuan, 
2007, 29, 989-994. 

[90] Mavaffi, M.D.; DiGirolamo, D.J.; Fan, Y.; Riddle, R.C.; Campbell, 
K.S.; van Groen, T.; Frank, S.J.; Sperling, M.A.; Esser, K.A.; Bam- 
man, M.M.; Clemens, T.L. Distinct growth hormone receptor 
signaling modes regulate skeletal muscle development and insulin 
sensitivity in mice. J. Clin. Invest., 2010, 120, 4007-4020. 

[91] Tamiya, T.; Kashiwagi, I.; Takahashi, R.; Yasukawa, H.; Yoshi- 
mura, A. Suppressors of cytokine signaling (SOCS) proteins and 
JAK/STAT pathways: regulation of T-cell inflammation by SOCS1 
and SOCS3. Arterioscler Thromb Vase Biol, 2011, 31, 980-985. 

[92] Rui, L.; Yuan, M.; Frantz, D.; Shoelson, S.; White, M.F. SOCS-1 
and SOCS-3 block insulin signaling by ubiquitin-mediated degra- 
dation of IRS1 and IRS2. J. Biol. Chem., 2002, 277, 42394-42398. 

[93] Cardinali, B.; Castellani, L.; Fasanaro, P.; Basso, A.; Alema, S.; 
Martelli, F.; Falcone, G. Microrna-221 and microrna-222 modulate 
differentiation and maturation of skeletal muscle cells. PLoS One, 
2009, 4, e7607. 

[94] Machova, P.K.; Lopotova, T.; Klamova, H.; Burda, P.; Trneny, M.; 
Stopka, T.; Moravcova, J. Expression patterns of microRNAs asso- 
ciated with CML phases and their disease related targets. Mol Can- 
cer, 2011, 10, 41. 



